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/~s , Abstract We review the Parisi-Frisch [T] MultiFractal formalism for Navier- 

Stokes turbulence with particular emphasis on the issue of statistical fluctu- 
ations of the dissipative scale. We do it for both Eulerian and Lagrangian 
Turbulence. We also show new results concerning the application of the for- 
malism to the case of Shell Models for turbulence. The latter case will allow 
us to discuss the issue of Reynolds number dependence and the role played 
by vorticity and vortex filaments in real turbulent flows. 
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ON i i Introduction 

Turbulent flows are characterized by strong fluctuations of energy transfer 
from small to large scales [2]. The best known theoretical approach to de- 
scribe the statistical properties of turbulent flows is due to Kolmogorov in 
f94f (K41). The K41 theory assumes that turbulent fluctuations can be 
considered homogeneous and isotropic at very small scales. Moreover, Kol- 
mogorov assumed that the average (in time) rate of energy dissipation e is 
Reynolds independent for large enough Reynolds number, Re = UqLq/v, 
where Uq is some characteristic velocity at large scale Lq and v is the kine- 
matic viscosity of the flow. For the sake of simplicity we will set Uq = Lq = 1 
from now on. Thus Re — v^ 1 . The equations of motions governing hydro- 
dynamical incompressible turbulence are the Navier-Stokes equations, with 
density set to p = 1 : 

d t u + (ud r )u = d r p + vAvl, (f) 
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where Ui is a three-dimensional field satisfying diUi — 0. Note that in terms of 
Ui the average rate of energy dissipation is given by e = v (pi u j) • Thus, 
e becomes independent of Re if at large Reynolds number diUj ~ Re 1 / 2 . 

The statistical properties of the vector field u, can be studied by using the 
n-order tensor, given by the correlation functions (iij 1 (x ai )...'Uj n (x 0ra )), where 
the index i\, ..i n stands for the spatial directions and x Qi are n-points in the 
physical space, and the symbol (...) means ensemble average. If we assume, 
following Kolmogorov, that turbulence becomes isotropic and homogeneous 
at scales small enough, then we can rewrite the correlation functions in terms 
of only longitudinal velocity increments, S r u — 5u(r) ■ f (where <5u(r) = 
u(x+r) — u(x)) and transverse increments, 5 r w = <5u(rx), (where r^-u = 0). 
For a discussion on the complexities arising when also anisotropic fluctuations 
are relevant see [3]. 

One of the main predictions of the Kolmogorov theory is contained in his 
famous 4/5-equation for homogeneous and isotropic turbulence: 

((S r u) 3 ) = -ier + 6is^-((5 r u) 2 ). (2) 

5 ar 

For large Re (i.e. small v), and at fixed distance r, the 4/5 equation tells us 
two important information: the statistical properties of S r u have a non zero 
skewness, and they have a scaling behavior as a functions of r. Note that 
(|2|) , taken in the limit of vanishing viscosity, is dimensionally consistent with 
having S r u ~ r 1 / 3 , which implies for the p-th order longitudinal Structure 
Functions, S (p \r) = {(S r u) p ) ~ r p / 3 . 

It turns out experimentally that the K41 theory is only partially correct, 
i.e. equation ([2]) is very well verified when the reference scale, r, is much 
smaller than the integral scale, L$, implying a recovery of homogeneity and 
isotropy for small enough scales. However, the structure function data are 
not in agreement with the Kolmogorov theory for all orders p. For instance, 
K41 theory predicts that the hyper-Flatness, S^ p \r) / '{S^ 2 \r)) 2 , should be 
constant as a function of r, while it is not. Both experiments and numerics 
show that it grows when going from large to small scales. This effect is usually 
referred to as "intermittency" and it is clearly missing in the K41 scenario 
(for a collection of experimental and numerical data on both Eulerian and 
Lagrangian data, see [HE]). 

What is the physical explanation of intermittency? Is intermittency a fi- 
nite Reynolds effect, i.e. does it disappear asymptotically for large Rel Are 
scaling properties of the turbulent field spoiled by intermittency? These ques- 
tions are the fundamental issues discussed by the scientific community in the 
last 50 years. Unfortunately, there are no exact results about intermittency 
in Navier-Stokes turbulence. So, we must rely both on experimental and nu- 
merical data, or on phenomenological and dynamical modeling of turbulence. 
In particular, a key problem is to control the effects of viscosity, i.e. how to 
connect experimental and numerical data - naturally limited in Reynolds 
numbers - to any theoretical or phenomenological understanding proposed 
for the asymptotic infinite Reynolds number case. In this paper we will deal 
with such an issue, focusing on the MultiFractal (MF) phenomenological de- 



3 



scription of "infinite Reynolds" turbulence and on its extension to the case 
of finite Reynolds. 

The paper is organized as follows. First, we review some of the results so 
far obtained by employing the MultiFractal theory of turbulence, introduced 
25 years ago by Parisi and Frisch [IJ[2] (see [6] for a recent review). Then, we 
discuss subtleties connected on how to introduce in a self-consistent way the 
effects of viscosity in the MF description, i.e. how to control finite Reynolds 
effects for both Eulerian and Lagrangian turbulence. Finally, we present new 
data obtained on a class of dynamical systems for the turbulent energy cas- 
cade, the so-called Shell Models [SHU] . This will allow us to address and test 
the MF formalism at changing Reynolds numbers. We conclude with some 
perspectives for further work. 

2 Review of the multifractal approach of turbulence 

The multifractal approach to turbulence is based on the assumption that the 
statistical properties of turbulent flows do exhibit scaling properties even if 
there is intermittency. Let us discuss scaling in a rather more abstract way by 
using the symmetries underlying Euler and Navier-Stokes equations. The Eu- 
ler equations have the "standard" symmetries of the Newton laws for inviscid 
fluid, i.e. isotropy, parity and time reversal. Moreover, the Euler equations 
display another important set of symmetries of global scale invariance, i.e. 
invariance under the transformations [9]: 

r^Xr; u -> A ,1 m; t -» X^H. (3) 

The value of h is arbitrary and it is not fixed by any physical and/or mathe- 
matical constraint. The Navier-Stokes equations have no time reversal (due 
to viscosity) and satisfy a new scale invariance property, namely: 

r^Xr: u -> X h u; t -» A 1_fc t; v A 1+ V (4) 

Note that the new scaling property implies e — > A 3,l_1 e. We can now rephrase 
the K41 theory in the following way: using |(5J) the energy dissipation rate 
should be constant and equal to S^(r)/r which implies h = 1/3. In other 
words, the scaling transformation r — > Xr (A <C 1) is equivalent to a change 
in the viscous effect and, therefore, to a change of e. If we assume that e 
is homogeneous in space and time, it follows that h cannot be arbitrary 
and should be fixed to h = 1/3. Let us note that the scale invariance for the 
Navier-Stokes equations is now restricted to a particular value of h, i.e. there 
is less "scale" symmetry in the system as naively predicted by the equation 
of motions. However, one can take a different point of view, namely one can 
assume that h is still arbitrary, i.e. the scale symmetry holds, although e is 
constant only on average. The problem is to understand what is meant by 
"average" . This is the conceptual step performed in the Parisi-Frisch paper 
PQ . The idea is to assume that in the limit of infinite Reynolds numbers, and 
for any fixed scale, r, the scale invariance holds with some probability P r (h). 
It is then assumed that 6 r u and Pft(r) are scaling function of r, i.e. 



6 r u~ r h ; P h (r)~r F W. 



(5) 
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Then, all the correlation functions should be computed by averaging over the 
probability Ph(r): 

S (p) {r) oc J dh r F ^ +ph . (6) 

The breaking of time reversal still leads to e ~ const which implies the 
constrain S^(r)/r ~ const. Using the definition: 

£(p) = mhi[F(h) + ph], 

h 

we can perform the integral by the saddle point method and obtain (r) ~ 
/•C(p) with the constraint £(3) = 1. 

Historically, F(h) has been written as 3 — D(h) by assuming that D(h) is 
the fractal dimension of the set at scale r where 8 r u ~ r h and, for this reason, 
the above picture has been named the "multifractal" approach to turbulence. 
We note that there is no reason to introduce any fractal dimension, i.e. there 
is no a priori geometrical interpretation in defining Ph(r). 

The crucial physical point, which is discussed in this paper, is the role 
of the viscosity in the multifractal approach. The scaling invariance of the 
Navier-Stokes equations tells us that by a scale transformation r — > Xr noth- 
ing changes if the viscosity becomes smaller by a factor A 1+ ' 1 . In the real 
world, the viscosity is fixed and therefore it increases as 

A-U+''). Eventually, 

the viscosity becomes so large as to kill any turbulent fluctuations. We can 
assume that viscosity is relevant if the turbulent characteristic time scale 
r/8v{r) is of the same order of the viscous time scale r 2 /v, which defines a 
fluctuating, i.e. /i-dependent, dissipation scale T)(h): 

^ ~ 1 -> n{h) ~ yVtt+A). (7) 

The idea of a fluctuating dissipation scaling has been introduced by Paladin 
and Vulpiani [TUj and it should be considered as a consequence of the scale 
invariance (|2]). Let us discuss in more details this point. 

A possible interpretation of the scale invariance would suggest that the 
velocity difference S r u - and this has implications for the n-point correlation 
functions - depends on h in the following manner: 

S r u ~ 9(r/vku)f(r/Vkii)) h , (8) 

where r/kn is the Kolmogorov -not fluctuating- dissipative scale estimated 
as r?fc4i = (^ 3 /e) 1 ^ 4 . The appropriate asymptotic behaviors for the functions 
g(x) and f(x) must be: g{x) ~ const and f(x) ~ x for x 3> 1, while g{x) ~ x 
and f(x) ~ const for i<1. Thus f{x) and g{x) can be interpreted as cutoff 
functions on the inertial-range scaling behavior introduced by the dissipation. 
Let us remark that we can also expect a cutoff function entering into the 
definition of the probability distribution Ph(r), i.e. we should expect that 

Ph{r) ~ [f(r/vk4i)f (h) . (9) 
It is interesting to remark that in the range of scales where 



g(x) ~ const. 



(10) 
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we have, using © that S^(r) ~ [r f(r/r))]^ p \ Recalling that C(3) = 1, the 
previous expression can be rewritten as S^(r) ~ (S (3 \r))^P\ It turn out 
that the range of scales where the equality (fT0|) holds is much more extended 
then the usual inertial range, i.e. from experimental and numerical data one 
can safely deduce that the relation (| 10(1 is valid down to scales x ~ 10; i.e. 
where viscosity is already important. As a result, both function f(x) and g{x) 
are not constant anymore at those scales, but the viscous corrections on both 
functions are almost the same. This empirical fact is known as ESS (Extended 
Self Similarity), it has been extensively used in literature [TT|rT2] to estimate 
the scaling exponents ((p) with good accuracy. Eq. (|10p cannot hold also for 
very small scales, i<Cl, where viscosity starts to have a dominant role. In 
order to take into account also of these effects, we need to consider relation 
(|7|) and its consequences. The generalization can be done following [13] 

S r u ~ g(r/v(hW h (r/r,(h)) h ; P h (r) ~ [f(r/v(h))f W . (11) 

In section (|3.2l) . we will discuss some explicit form of the function f(x) and 
g(x). Equations (fTTjl tells us how to introduce in a self-consistent way the 
effect of finite viscosity, i.e. breaking of scale invariance) within the scale in- 
variant MF approach. Wc will refer to the whole picture as the multifractal 
conjecture. More precisely, we must expect that for all quantities that are 
invariant with respect to the Navier-Stokes symmetries, the scaling proper- 
ties can be obtained by using the same universal function F(h), including 
the effect of finite viscosity. Thus, all the Reynolds number effects on the 
statistical properties of turbulence should be predicted using the function 
F(h), the scaling invariance and equation ([7|). From this point of view, the 
multifractal conjecture should be able to predict how different quantities 
(for instance moments of the velocity gradients) depend on Re. The explicit 
functional form for the interpolating functions f(x) and g{x) may be safely 
considered not particularly relevant, as soon as the asymptotic scaling limits 
are preserved (see also [l4l[T5] for some preliminary attempts to determine 
them using constraints from the Navier-Stokes equations). 

Let us make some remarks on the above picture, mostly in a historical 
perspective. From the theoretical point of view, the major challenge is to 
compute the function F{h). In some particular, linear cases, this computation 
can be done [16] . although no one knows how to perform such calculations 
for the Navier-Stokes equations. There has been a lot of effort to clarify, 
using both experiments and numerical simulations, whether the statistical 
properties of turbulence exhibit scaling behavior. So far, there is a wide 
consensus that this is the case. Knowing the exponents £(p) one can easily 
compute the function F(h) by some reasonable "fitting" procedure. Thus, we 
can verify whether the multifractal conjecture is satisfied or not, i.e. whether 
there exists or not some quantity (invariant with respect the group of Navier- 
Stokes symmetry) which cannot be predicted knowing the function F(h). For 
non-invariant quantities, the multifractal conjecture tells us nothing. Also, 
there have been some approaches in the past claiming that for large enough 
Re, the statistical properties of turbulence are asymptotic to the K41 theory. 
Experiments and numerical simulations do not support such claims [4ll5|ll7j , 
so far. 
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3 Predictions concerning the fluctuating dissipative scale 

3.1 Eulerian Framework 

The multifractal conjecture allows us to develop a number of different pre- 
dictions. Let us assume that we know the function F(h) or, equivalently, the 
scaling exponents C(p)- It is then possible to predict the scaling behavior of 
the following non trivial quantities. Some of these predictions are based on 
the basic idea that the dissipation scale of turbulence is a fluctuating quantity 
satisfying eq. (jTj) . It is then possible to compute the probability distribution 
of the velocity gradients by using the estimate [1011181119] : 

o r u ~ . (12) 

V 

Since 8 n (h) u ~ { J l{h)) h and rj(h) ~ v 1 /( 1+h \ we have: 

((d r u)") ~ I dhvMh-y+FWmi+h) „ifex« (13) 



which predicts the scaling properties of the gradient as function of Re ~ 
and where the last expression is obtained by the saddle point method when 
the Reynolds number Re ~ v~ l tends to infinity. Here 

x (q) = min[(2g (h - 1) + F(h))/(1 + h)}. 

h 

Note that the prediction is highly non trivial and it is connects x(p) with 
the exponents £(p) via a non-linear relation valid for any function F(h), (see 

my- 

, \ p-C(p) .,, C(p)+p nA , 
x(q) = — 2 — ; Wlth q = — 2 — ' *■ ' 

The above relation implies the existence of dissipative anomaly, i.e. x(2) = 1, 
for any spectrum such that £(3) = 1: 

lim e ~ Re" 1 ((d r u) 2 } — > const. (15) 

Re — >oc 

Such non-trivial multifractal spectrum for velocity gradients may be used to 
define a multifractal spectrum of dissipative scales, each-one defining the dif- 
ferent cross-over from the inertial range to the dissipative range for different 
moments. To implement this, we may follow two slightly different methods. 
First, let us define the n-th order dissipative scale as the intersection be- 
tween the extrapolation for large scales of the differential smooth behavior 
of S^ n \r) with the extrapolation for small scales of the inertial behavior: 

SW(r)~r"((cU) n ) if r« 77,41 (16) 



S(»)( r ) ~r c(jl) if r»Tfc4i. 



(17) 
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The intersection between the two power-law behavior can be considered a 
good estimate for the typical dissipative scale, r)( n -) , of the n-th order structure 
functions: 

( V(n) T((d r uT) = fa (n) )« W >. (18) 

Using expression ([14"]) for the scaling of gradients, we get: 

77 (n) ~ i? e -x(«)/(«-C(»)). ( 19 ) 

Notice that according to this definition we have for the dissipative cut-off of 
the second order structure function: 

77(2) ~ ife-x(3)/<2-C(3» (20) 

Another, slightly different, way to define the moment-dependent viscous cut- 
off is the one followed in refs. pDl[2Tj . where a balancing inside the equations 
of motion between inertial and dissipative terms is used to obtain r)r n -\ . Fol- 
lowing [20,21J, it is easy to derive the dimensional balancing between the 
inertial terms and the dissipative contribution in the evolution of the n-th 
order generic structure functions, from (TT|): 

S (n+1) (r)/r ~ vS {n) {r)/r 2 . (21) 

If the n-th order dissipative scale is now obtained by asking it to be the scale 
where the previous matching holds, we easily obtain, another expression for 
the cut-offs [20] : 

r)[ n) ~i? e -V(C(»+D-C(n)+i). (22) 

Let us notice that the two expressions ([19} and (|22| are exactly equivalent for 
the cut-off entering in the dissipative anomaly 77(2) = rjL-, = Re^ 1 ^ 2 ^ l '^\ 
for any choice of F(h), while they are slightly different for high-order mo- 
ments. The discrepancies are however very small, as can be directly checked 
by plugging some possible shape of the £(p) exponents that fit the experi- 
mental data (for example the log-Poisson expression derived in [22] , as shown 
in fig. [J). Recently, ad- hoc advanced computations have been performed by 
changing the numerical resolution, in order to test the existence of such non- 
trivial fluctuations at dissipative scales [23]. The numerical data have also 
been succesfully compared with the MF prediction ([7]) in [24]. Note that 
it is quite difficult to get good laboratory data for velocity fluctuations at 
scales much smaller than 77^41 due to the intrusive nature of the experimental 
probes. Nevertheless, in recent years very interesting and promising experi- 
mental techniques have been developed to track particles in turbulent flows, 
accessing temporal fluctuations instead of spatial fluctuations |25lf26l[27l[28l 
I3UII29"] . The quality of predictions based on the multifractal conjecture is even 
more striking by performing the statistical analysis of the velocity field in the 
Lagrangian framework. 
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3.2 Lagrangian Framework 

Given the importance of particle dynamics in turbulent flows, numerous nu- 
merical and experimental studies have flourished in the last few years |26[ 
I3DI32. 33, 34, 35, 36 (see also [3D] for a recent review). Neutrally buoyant par- 
ticles, advected by a turbulent velocity field u(x, t), follow the same path 
as fluid molecules and evolve according to the dynamics X(t) = v(t) = 
u(X(i),i), where the Lagrangian velocity v equals the Eulerian one u com- 
puted at the particle position X. Such particles constitute a clear-cut indi- 
cator of the underlying turbulent fluctuations. Recently, it has been shown, 
by comparing the different numerical studies and different experimental re- 
sults 5,37 , that Lagrangian turbulence is universal, intermittent, and well 
described by a suitable generalization of the Eulerian Multifractal formalism 
to the Lagrangian domain [38 , 39 , 40 . 

Lagrangian Structure Functions (LSF) are defined as: 

5g(r) = (Mt + r) -«,(<)]*) = ((8 T ViY), (23) 

where i — x,y, z runs over the three velocity components, and the average is 
defined over the ensemble of particle trajectories evolving in the flow. From 
now on, we will assume isotropy and therefore drop the dependency from the 
spatial direction i. The presence of long spatial and temporal correlations 
suggests, in analogy with critical phenomena, the existence of scaling laws 
for time scales larger than the dissipative Kolmogorov time and smaller than 
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the typical large-scale time, t v <C t <C Tl:xs 

5^ ) (r)~T* (p) . (24) 

Straightforward dimensional arguments d Za Kolmogorov predict z(p) = p/2, 
independently of the flow properties. However, it is known that LSF expe- 
rience show strong variations when changing the time lags r, as highlighted 
by the increasingly non-Gaussian tails characterizing the probability density 
functions of S T v for smaller and smaller r's [26 . This leads to a breakdown 
of the dimensional argument; correspondingly there is a growth of the La- 
grangian flatness when going to smaller and smaller time lags In [51137) it 
has been shown that the scaling exponent re(4) = z(4)/z(2), entering in the 
evolution of the fourth-order flatness: 



(4 2) (r)) 



2 



((S ( 2 ) (r)r (4)-2 (25) 



does not depend on the experimental or numerical large-scale set-up. In other 
words, the high frequency fluctuations are universal. 

It is possible to obtain a link between Eulerian and Lagrangian MF for- 
malism via the dimensional relation |38U39j : 

t ~ r/S r u. (26) 

Indeed, if we substitute (|26|) into (O and into ([6]) we obtain a Lagrangian 
prediction for LSF once the Eulerian one is known via its F(h) spectrum 
[39], namely 

4 P) (t) ~ J dhT^ h+F ^'^ ~ T*<*> , (27) 

where 

z{ P ) = mm[{ph + F{h))/{1 + h)}. 

h 

This result is obtained by the saddle point method for inertial-range time 
intervals: t v < t < Ti. Let us notice that the above relation (j2"T|) can be 
read as a prediction for the Lagrangian domain, once the Eulerian statistics 
is known. This is because we are using the same MF functions F(h) for 
both domains. The suggested road map should be the following: (i) first 
measure the Eulerian scaling exponents C(p)j (ii) then via an inverse Legendrc 
transform extract the F(/i)-spectrum; (iii) finally, apply the relation (f27|) 
and calculate the Lagrangian scaling. Such procedure is working well, at 
least within the statistical limitation and the Reynolds number limitations 
allowed by numerical and experimental state-of-the-art techniques [S] (see 
also [H] for a recent theoretical attempt). 

Even more interesting, in [5], the same argument, leading to the spa- 
tial dissipative fluctuating scale 0, has been extended into the Lagrangian 
domain to obtain an expression for the fluctuating dissipative time scale [35] : 

r v {h) ~ „(i-'0/(i+'0 . (28) 
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Following the same ideas discussed in Sec. [2] about the viscous modification 
to the MF formalism, we may now introduce two functions, f{^jj^) and 
.9(7^)), which take into account the viscous corrections to the Lagrangian 
inertial scale (|27|l. Is is customary to use for the two cross-over functions 
a Batchelor-Meneveau functional form [13,,4Ql[5]. The global description for 
velocity increments becomes then: 

S T v(h, T v (h)) T -l^- — , (29) 

in which [3 is a free parameter which controls the crossover from dissipative 
to inertial time lags. It is easy to realize that the above expression reproduces 
the asymptotic regimes, 5 T v ~ r for r <C t v and S T v ~ t 11 ^ 1 ^^ in the inertial 
range. As in the Eulerian case, each exponent h must be weighted with its 
probability, which is also modeled to include dissipative-range physics, mim- 
icking what has already been discussed in the introduction for the Eulerian 



P h (T,T 11 (h))=Z-\T) 



(30) 



Here Z is a normalizing factor. We notice that again everything is expressed 
in terms of the same MF function F(h). Only a new free parameter, /3, 
enters to describe the whole Lagrangian behavior, at all time scales, once 
the Eulerian MF spectrum, F(h), is known. Putting everything together, we 
obtain for the LSF the following MF prediction: 



((S T v)P) ~ / dhP h (T,T v )[5 T v(h)F . (31) 



The previous formula for the Lagrangian domain predicts a highly non trivial 
shape for the local scaling exponents entering in the Flatness behavior: 



d log(4 2 >(r)) 

As shown in fig. ([2]) , where we plot it for a series of different Reynolds num- 
bers and for a given choice of the multifractal spectrum F{h) (see caption 
for details). The functional dependence shown in fig. @ was tested in [S] 
only for k(4) and for a limited set of Reynolds numbers, due to the unavoid- 
able limitations in experiments and numerical simulations. We want now to 
test it further, by boosting the Reynolds number by orders of magnitudes. 
To achieve this, we switch to a dynamical surrogate of the Navier-Stokes 
equations. 
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Fig. 2 Fourth-order local scaling exponents, k(4, t) 

formalism |g7J for three different Reynolds numbers, 7?e = 2.10 4 , 9.10 4 , 7.10 5 (from 
top to bottom) The two horizontal line correspond respectively to the non intermit- 
tent, K41 case, «(4) = 2, and to the infinite Reynolds number limits, k(4) = 1.72. 
The F(h) spectrum used to derive this values has a log-Poissou shape as proposed 
in [22]. The value of the free parameter is /3 = 4. Notice that the inertial range 
extension, identified as the region where the local exponent is constant, becomes 
larger and larger when increasing the Reynolds number. Notice also the strong in- 
crease in the intermittency, measured by the deviation from the K41 value, k — 2, 
for time scales across the viscous domain. 



3.3 Shell Models 



Let us now discuss in more details the behavior of the local scaling exponents 
for the Lagrangian structure functions (|3^|) . A decrease of the local scaling ex- 
ponent simply means that intermittency is growing in the dissipation range. A 
detailed analysis of the flow configurations, which can be performed in high- 
resolution numerical simulations, show that the tip in the scaling exponents 
is strongly related to the presence of coherent three-dimensional vortices [34j 
143] . This is an important observation which needs to be investigated care- 
fully. The physical question we are discussing, concerns the possibility that 
coherent structures, i.e. vortices, are responsible for intermittency in three- 
dimensional turbulence, not just in the dissipative region but for the whole 
inertial range. After all, if the largest intermittent fluctuations of velocity 
gradients are due to vortices, it is somehow reasonable to think of the effect 
of coherent structure as a key feature for explaining intermittency. If this is 
true, the multifractal conjecture may be misleading, i.e. scale invariance is 
not a fundamental properties of three-dimensional turbulence. Nevertheless, 
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in [5j it has been shown that the multifractal conjecture does predict the tip 
in the local scaling exponents for the structure functions. 

So far, we have discussed the statistical properties of homogeneous and 
isotropic turbulence by checking the multifractal predictions against numer- 
ical and experimental data. In the last twenty years, it has been shown that 
there exists a class of simplified models, named shell model, which shows mul- 
tifractal intermittency (anomalous scaling) similar qualitatively and quanti- 
tatively to what it is observed in Navier-Stokes turbulence. Among many 
different shell models, we shall consider the shell model proposed in [H] ( see 
also [8] for a review). Shell models of turbulence, can be seen as a truncated 
description of the Navier-Stokes dynamics, preserving some of the structure 
and conservation laws of the original equations but destroying all spatial 
structures. They are described by the following set of ODEs: 



i(k n+ iu* l+1 u n+2 - 5fc„u*_ x u n+ i 

+ (1 - <5)fc n _i1i n _iU n _ 2 ) + fn ■ 



(33) 



Here the u n s are the velocity modes restricted to 'wavevectors' k n = fco2" 
with k determined by the inverse outer scale of turbulence. The model con- 
tains one free parameter, 5, and it conserves two quadratic invariants (when 
the force and the dissipation term are absent) for all values of S. The first 
is the total energy J^ ra \u n \ 2 and the second is a sort of generalized helicity 
Sn( — \ u n\ 2 , where a = log 2 (l — S) [8]. Here we consider values of the 
parameters such that < S < 1. The scaling exponents characterize the shell 
model structure functions, defined as 

#W = K>~^ (2) , (34) 
(k n ) = Im (u n ^u nU ; i+1 ) ~ fc-^ 3 ' , (35) 
Sp{k n ) ~ k n ^ . 

The values of the scaling exponents have been determined accurately by di- 
rect numerical simulations. Besides £(3) which is exactly unity, because a 
relation similar to ^ holds, all the other exponents £(p) are anomalous, 
differing from p/3. For 8 = —0.4, the value of the ((p) are close to the scal- 
ing exponents of the Navier-Stokes equations. In applying the multifractal 
conjecture to the shell model, we shall assume u n <~ k~ h with probability 
kn F< " h \ while the fluctuating dissipative scale kd(h) is defined by the relation 
Ud/{kdv) ~ 1, where — k^ h . In the dissipation range, the behavior of the 
shell model is roughly consistent with u n ~ k n exp(— k n /kd)^ By matching 
between the inertial range (i.e. u n ~ k~ h ) and the dissipation range, we then 
obtain: 

u n ~ k~ h (l + A(k n /k d (h)) 1+h ) exp(-k n /k d (h)) (36) 

1 Balancing of the nonlinear and viscous terms in the far dissipation range 
actually gives to leading order u n ~ fc n exp( — k") with a ~ 0.69 solution of 
2 "« + 4" a = 1. 
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where A is a Reynolds independent quantity. Consequently, the probability 
distribution Ph(k n ) should be modified in order to take into account the 
dissipation effects, namely: 

Ph(k n ) ~ Ml + A{k n /k d {h)))- l ]- F ^. (37) 

Knowing F(h), which as usual can be estimated from the knowledge of C(p)i 
we can compute the k n dependency of S^(k n ) both in the inertial and in 
the dissipative range and compare our findings with numerical simulations 
of P|l. 

In fig. ([3]) we plot the local scaling exponents 

K(4 ,M- dlog(Mkn) 



d io g (s 2 (M) 



computed from the numerical simulations of the shell model (upper left panel) 
and predicted by the multifractal conjecture (upper right panel). Different 
symbols refer to different Reynolds numbers. The first striking result is that 
the numerical simulations of the model clearly show a well defined tip in the 
dissipative region (i.e. large value of n in the figure), similar to what it is 
observed in experiments and in the numerical simulations. The tip is increas- 
ing towards small scales as the Reynolds number increases and it deepens. 
The straight line in the figure is a qualitative fit on the behavior of the tip 
as a function of Re. Clearly, the increase of intcrmittency in the dissipative 
region is scaling as log(i?e). Notice that this scaling behaviour was not vis- 
ible in the experimental and numerical data shown in [51137] because of the 
limited range of Reynolds spanned in those cases. In the right panel, we show 
k(4, k n ) as predicted by using the multifractal conjecture with A — 0.8. As 
one can observe, the qualitative and quantitative predictions are in remark- 
able agreement with the numerical simulations. 

The quality of the result does not change by increasing the order of the 
structure functions. In fig. we also show behavior of k(6, k n ) (i.e. the 
local scaling exponents of the sixth order structure functions) and compare 
the numerical simulations with the multifractal prediction. Again, the results 
confirms what we found for the lower order. Notice again, the widening of 
the inertial range when increasing the Reynolds number and the power-law 
behavior of the wavenumber where we observe the highest intcrmittency, i.e. 
the minimum in the local scaling exponents. 

From the above analysis we can draw some interesting conclusions. First 
of all, the increase of intermittency in the dissipative range is not due to 
coherent vortices (there are no vortices in the shell model). Moreover, the 
increase of intermittency is predicted by the multifractal conjecture because 
of the fluctuations of the dissipative scale. Translating back these results to 
the real Navier-Stokes equation, we are tempted to conclude that coherent 
structures exist but their dynamics is nor relevant to explain intermittency 
in turbulent flows. They are the tail more than the dog (see chapter 8.9 of 
[2]). Actually, one can take the opposite point of view: the matching between 
inertial range (scale invariance) intermittency and the dissipation range pro- 
duces an increase in fluctuations and, consequently, an increase of vorticity. 
This effect is dominant at low Reynolds number, as it is clearly observed in 
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5 10 15 20 5 10 15 20 

l°92(k n ) >°92(K) 



Fig. 3 Top: local scaling exponents for fourth order structure functions, k(4, k n ), 
in the Shell Model, versus log 2 (fc n ). Left numerical results at four Reynolds num- 
bers in the range Re £ [10 2 : 10 4 ]. Right, the MF prediction using relations (|36l37Jl ; 
the straight line is a fit for the behavior of the bottleneck (maximum of intermit- 
tency), the slope of the line is —0.028. Bottom: the same for the sixth order scaling 
exponents k(6, k n ). The straight line has the slope —0.056 

the shell model: i.e. at low Reynolds numbers no scaling behavior is observed 
and intermittency is strongly dominated by the fluctuations in the dissipation 
range. Let us also point out that our observations is in qualitative agreement 
to the known phenomenology of boundary layer turbulence. Near the wall, 
the local Reynolds number is relatively low and strong intermittency is ob- 
served together with a rich dynamic behavior of coherent structures (hairpin 
vortices) . The full dynamics is dominated by strong intermittent fluctuations 
and it is tempting to relate our previous discussion to this specific well known 
turbulent flows. 

What is remarkable in the above discussion, is that the change in the 
intermittency level, is explained by the same MF spectrum F(h), at all scales 
and at all Reynolds numbers. In other words, the study of low Reynolds 
turbulence, and of the transition from viscous to inertial range 45,46J, may 
teach us a lot with respect to the high Reynolds limits [47] . 



4 Conclusions and Perspectives 

Let us summarize our main points. First, we have reviewed the ideas about 
how to introduce, in a self-consistent way, dissipative effects in the Multi- 
fractal description of turbulence, both Eulerian and Lagrangian. The MF 
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formalism, predicts an enhancement of intermittency in the so-called inter- 
mediate viscous range [48], as measured by the local scaling exponents (see 
fig, [2]). We have commented on the fact that such a trend is in very good 
agreement with real Lagrangian Turbulence data [3] , at least concerning low 
order moments and moderate Reynolds numbers, i.e. up to the numerical and 
experimental state-of-the-art. Second, in order to test the formalism also for 
high Reynolds numbers, we switched to a class of Shell Models of turbulence. 
Here, thanks to the much simpler structures of the model, the Reynolds de- 
pendence of the MF prediction can be also tested. We have specialized our 
discussion on the enhancement of intermittency measured around dissipative 
scales. The existence of the very same phenomenon observed in real Navier- 
Stokes eqs, lead us to conclude that the viscous intermittency is not due to 
vortex filaments. Many problems remain opened. For example, the Batchelor- 
Meneveau structure presented in (|29j) is not compatible with the requirement 
that velocity fluctuations follow a simple multiplicative local -in scales- pro- 
cess from large to small scales. This is due to the fact that the fluctuating 
temporal scale T v (h) appears in the definition of the velocity increments and 
it depends on them (in eqs. (|28ti29[l the local scaling exponent h is the same). 
In other words, the functional relation (|29|) introduces non-local correlation 
between inertial and dissipative scales. A local-in-scale multiplicative process 
which takes into account the fluctuating cutoff can be introduced somewhat 
empirically by building a multiplicative cascade and stopping it according to 
the criteria J7J . It would be interesting to see if such a procedure reproduces 
the Navier-Stokes data and the Shell Model data as nicely as (|2U1) . 

We acknowledge many fruitful collaborations on these issues with M. 
Cencini, A.S. Lanotte and F. Toschi. L.B. wants to acknowledge also fruitful 
discussions with V. Yakhot. 



References 

1. G. Parisi and U. Frisch in Turbulence and Predictability in Geophysical Fluid 
Dynamics, Proceed. Intern. School, of Physics "E. Fermi", eds. Ghil M, Benzi 
R, Parisi G. 84-87. (1985). 

2. U. Frisch, Turbulence: the Legacy of A.N. Kolmogorov (Cambridge University 
Press, Cambridge UK, 1995). 

3. L. Biferale and I. Procaccia, Phys. Rep. 414, 43 (2005). 

4. A. Arneodo et al., Europhys. Lett. 34, 411 (1996). 

5. ICTR Coll.; A. Arneodo et al., Phys. Rev. Lett 100, 254504 (2008). 

6. G. Boffetta, A. Mazzino and A. Vulpiani. J. Phys. A 41 363001 (2008). 

7. T. Bohr, MH Jensen, G. Paladin and A. Vulpiani. Dynamical systems approach 
to turbulence. (Cambridge University Press, Cambridge, UK, 1998) 

8. L. Biferale Ann. Rev. Fluid Mech. 35, 441, (2003). 

9. U. Frisch in Turbulence and Predictability in Geophysical Fluid Dynamics, Pro- 
ceed. Intern. School, of Physics "E. Fermi", eds. Ghil M, Benzi R, Parisi G. 71 
(1985). 

10. G. Paladin and A. Vulpiani, Phys. Rev. A 35, R1971 (1987). 

11. R. Benzi et al., Phys. Rev. E 48, R29 (1993). 

12. R. Benzi et al. Physica D 96, 162, (1996). 

13. C. Meneveau, Phys. Rev. E 54, 3657 (1996). 

14. L. Chevillard, B. Castaing, E. Leveque and A. Arneodo. Physica D 218, 77 
(2006). 

15. L. Sirovich, L. Smith and V. Yakhot. Phys. Rev. Lett. 72, 344 (1994). 



16 



16. K. Gawedzki and A. Kupiainen. Phys. Rev. Lett. 75 3834 (1995). 

17. T. Gotoh, D. Fukayama and T. Nakano, Phys. Fluids 14, 1065 (2002). 

18. M. Nelkin, Phys. Rev. A 42, 7226 (1990). 

19. R. Benzi, L. Biferale, G. Paladin, A. Vulpiani and M. Vergassola. Phys. Rev. 
Lett. 67, 2299, (1991). 

20. V. Yakhot and K. R. Sreenivasan. Physica A 343, 147 (2004). 

21. V. Yakhot. Physica D 215, 166 (2006). 

22. Z-S. She and E. Leveque, Phys. Rev. Lett. 72, 336 (1994). 

23. J. Schumacher. Europhys. Lett. 80, 54001 (2007). 

24. L. Biferale, Phys. Fluids 20 , 031703 (2008) 

25. A. La Porta et al., Nature 409, 1017 (2001). 

26. N. Mordant, P. Metz, O. Michel and J.-F. Pinton, Phys. Rev. Lett. 87, 214501 
(2001). 

27. S. Ott and J. Mann, New J. Phys. 7, 142 (2005). 

28. B. Luthi, A. Tsinober and W. Kinzelbach, J. Fluid Mech. 528, 87 (2005). 

29. S. Ayyalasomayajula, A. Gylfason, L.R. Collins, E. Bodenschatz and Z. 
Warhaft. Phys. Rev. Lett. 97 144507 (2006). 

30. F. Toschi and E. Bodenschatz, Ann. Rev. Fluid Mech. (to appear in 2009). 

31. J. Berg et al., Phys. Rev. E 74, 016304 (2006). 

32. H. Xu et al., Phys. Rev. Lett. 96, 024503 (2006). 

33. H. Homann et al. J. Plasma Physics 73, 821 (2007). 

34. L. Biferale et al., Phys. Fluids 17, 021701 (2005). 

35. P.K. Yeung, S.B. Pope and B.L. Sawford, J. Turb. 7 (58), 1 (2006). 

36. R. Fisher et al., IBM Journ. Res. Devel. 52 Num. (1/2) to appear (2008). 

37. L. Biferale, E. Bodenschatz, M. Cencini, A. Lanotte, N. Ouellette, F. Toschi 
and H. Xu. Phys . Fluids 20, 065103 (2008). 

38. M.S. Borgas, Phil. Trans. R. Soc. London A 342, 379 (1993). 

39. G. Boffetta, F. De Lillo and S. Musacchio, Phys. Rev. E 66, 066307 (2002). 

40. L. Chevillard et al., Phys. Rev. Lett. 91, 214502 (2003). 

41. K. P. Zybin, V. A. Sirota, A. S. Ilyin, and A.V. Gurevich. Phys. Rev. Lett. 
100, 174504 (2008). 

42. L. Biferale et al., Phys. Rev. Lett. 93, 064502 (2004). 

43. J. Bee et al. Phys. Fluids 18, 081702 (2006). 

44. V. L'vov, E. Podivilov, A. Pomyalov, I. Procaccia and D. Vandembroucq. Phys. 
Rev. E 58 1811 (1998). 

45. D. Lohse. Phys. Rev. Lett. 73 3223 (1994). 

46. D. Lohse and A. Mueller Groeling. Phys. Rev. Lett. 74 1747 (1995). 

47. J. Schumacher, K.R. Sreenivasan and V. Yakhot. New J. Phys. 9, 89 (2007). 

48. U. Frisch and M. Vergassola, Europhys. Lett. 14, 439 (1991). 



